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ABSTRACT 

A design is presented for on elec- 
trical analog computer which can 
solve non-steady state heat transfer 
problems in an extensive radial for- 
mation containing a moving cylin- 
drical source. The computer is used 
to simulate a simplified thermal oil 
recovery process in which heat trans- 
fer from a moving, constant tem- 
. perature source takes place radially 
by con d itci ion only. Temperature 
distribution curves are shown for 
several different assumed modes of 
travel for the heat source. The data 
are used to estimate the residual 
fuel requirements necessary to main- 
tain a self-propagating isothermal 
front for the particular system being 
studied. 

Although the computer is design- 
ed to represent an unique system, 
conversion factors can be adjusted 
to show the effects of changing the 
assumed values of system constants, 
such as thermal conductivity of the 
formation and temperature of the 
source isotherm. Examples are given 
to illustrate the relative influence the 



Ortsrimil mn.nu.script received in Petroleum 
Branch office on July 20, 1!>r;r>: revise* I manu- 
script received on Nov. 11, ID:! 5. Paper pre- 
sented at Petroleum Branch Fall Meetin* in 
New Orleans Oct. 2-3 and at Southern Cali- 
fornia Petroleum Station Fall Meeting in Los 
'Vnwles Oft. 20-21. 



vol.* 20i, i r 



assumed values of these variables 
will have upon the quantity of fuel 
required to maintain the source. 

The maximum effect of heat 
transfer by the flowing air stream 
is estimated mathematically, and 
shown to be markedly reduced. 

Inasmuch as the examples discuss- 
ed refer to a simplified process, the 
results are not directly applicable to 
a practical field operation. It is be- 
lieved, however, that the data illus- 
trate some general trends which are 
important hi thermal recovery proc- 
esses. 

Data from laboratory and field 
experimentation can be used to 
modify the computer to take into 
account the influence of heat trans- 
fer by various mechanisms, for ex- 
ample, by infected air and by fluids 
ahead of the hot zone. With such 
refinements it would be possible to 
estimate more accurately the limit- 
ing values for the rate of travel 
of the high temperature front and 
for the required air injection rates. 

Use of this computer should aid 
in evaluating the economic feasibil- 
ity of oil recovery by various 
thermal processes. It is hoped this 
paper will stimulate further work 
by others to help accomplish this 
objective. 



INTRODUCTION 

As the costs of exploration for 
and development of new oil fields 
rise, increasing the recovery of oil 
from established fields becomes more 
and more important. In fields con- 
taining low gravity crude oil, ulti- 
mate recoveries by primary means 
are sometimes less than 10 per cent 
and often total only 10 to "20 per 
cent of the oil in place. Conventional 
second a ry recovery methods are not 
always effective in increasing this 
total. Methods which have been 
proposed for increasing the recovery 
of viscous oils include thermal proc- 
esses for increasing the mobility of 
the oil. One method is to provide a 
moving heat source generated within 
the formation to heat the oil and a 
flowing gaseous medium (for ex- 
ample, air) to sweep the oil to "the 
producing wells. The possibility of 
such a process has intrigued engi- 
neers for many years, 1, = * 3 but little 
or no information has been available 
in the literature relative to the tech- 
nical or economic feasibility of this 
method. Recently, however, interest 
has been revived by publication of 
laboratory and field test results indi- 
cating the feasibility of maintaining 

3 References Riven at end of paper. 



a self propagating heat source in an 
oil formation. Very high ultimate 
recoveries were indicated to be pos- 
sible 4 - s -\ F 

To date, numerous laboratory 
tests have been devoted to studying 
the many variables involved, such 
as oil and oil-sand characteristics, 
pressure, and air flux. Of necessity, 
such tests are usually made in linear 
cells. The application of laboratory 
results to a radial system of suffi- 
cient scale to evaluate the economic 
feasibility of a thermal oil recovery 
process is extremely expensive and 
subject to considerable trial and 
error. One intermediate step is to use 
mathematical and analog methods to 
study means of utilizing the experi- 
mental information in the most ef- 
fective manner. 

To a large extent, the evaluation 
of economic feasibility depends upon 
a knowledge of such factors as sweep- 
out efficiency, well spacing, and air 
injection rates and pressures. A study 
of areal sweep of a combustion pat- 
tern has been reported in a recent 
paper by H. J. Ramey, Jr., and G. 
W. Nabor. T 

It is proposed in the present paper 
to describe a design for an electrical 
analog computer which can calculate 
the temperature profile resulting from 
non-steady state heat transfer from 
a moving source in an extensive ra- 
dial system. A method is presented 
for using such profiles to estimate 
the fuel required to maintain a mov- 
ing heat source. The effect of changes 
in the assumed values of different 
variables is demonstrated through 
changes in scale factors in the com- 
puter. 

In the present study the computer 
was used to simulate a simplified ther- 
mal process in which the single fac- 
tor of heat transfer by conduction 
alone was considered. Inasmuch as 
the nature of the thermal oil recov- 
ery process may predetermine the 
quantity of fuel available for generat- 
ing the heat source, and the heat 
transfer characteristics are dependent 
upon the physical operating condi- 
tions, it may be concluded that the 
optimum mode of travel of the source 
will be affected by additional factors 
not included in the present design. 
Through the utilization of data ob- 
tained in laboratory and field experi- 
ments, the computer described herein 
may be modified to evaluate different 
hypothetical mechanisms for the re- 
covery of oil by thermal means. 

APPARATUS AND PROCEDURE 

Mathematical or experimental so- 
lutions of unsteady state flow prob- 



lems are often difficult and time 
consuming. A common method of 
obtaining a rapid solution to such 
problems is often found by use of 
electrical analogs.' 810 One such net- 
work 11 has been adapted in our labo- 
ratory to study heat conduction from 
a moving, cylindrical heat source in 
the earth. 

The electrical model inakes use of 
the similarities which exist between 
the flow of heat in a rigid body and 
that of charge in a non-inductive elec- 
trical circuit. 

The method is based on the iden- 
tity in form between the following 
electrical and heat flow equations: . 
Heat Flow Electrical Flow. 
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It follows from this identity in 
form of the fundamental defining 
equations that all mathematical re- 
sults based on them will be similar. 

In the analog computer discussed 
in this paper, the simulated forma- 
tion was divided into 100 concentric, 
cylindrical tubes of unit thickness! 
each varying from the next one by 
an increment of radius. Because of 
the symmetry of this two-dimensional 
problem, it was possible to simplify 
the electrical network by adjusting 
the values of the components to rep- 
resent the thermal capacitance and 
resistance associated with the radial 
location of each thermal element. To 
approximate an infinitely large me- 
dium, the last element in the circuit 
was designed to correspond to a for- 
mation volume which was large com- 
pared to the entire radial system 
being studied. 

In order that the computer could 
be operated at a convenient voltage 
and transient time interval, scale fac- 
tors were chosen for the conversion 
from thermal to electrical units, as 
follows; 

/ = 0.18 v/°F; that is, 180 v is 

equivalent to 1000°F. 
n = 0.006 sec/hour; that is, one 
minute in the computer is 
equivalent to 10,000 hours 
actual time, 
m = 3.04 X 10"" farads/Btu/°F; 
that is, one microfarad is 
equivalent to 3290 Btu/°F. 
The fallowing thermal constants 
were assumed to be representative of 
a fluid saturated, porous rock: 

thermal conductivity, k = 1.6 

Btu/sq ft hour °F/ft 
specific heat, c = 0.2 Btu/lb °F 
density, p = 140 lb/cu ft 
diffusivity, « = 0.057 sq ft/hour 



The scale factors may easily be 
changed without invalidating the re- 
sults as long as the relationships be- 
tween the two systems (electrical and 
thermal) are maintained. It is possi- 
ble, within certain practical limits, to 
study the effect of changing the as- 
sumed values of the variables in the 
thermal systems by making suitable 
changes in the scale factors. This is 
important in considering the useful- 
ness of the computer. 

For a given position (x) in the 
thermal system, the thermal resist- 
ance is: 

Ai^v-ji 

= 0.625^ ° Fh ° Ur 
A x Btu 

To convert to electrical units: 
m 

L233 X 10' ^£ ohms 

where A* is the logarithmic mean 
area. 

For the same element, the thermal 
capacitance is: 

AC t = 7r(r* +i - r) - 1 • P • c 
= 87.6 0* + - r) Btu/«F 

To convert to electrical units: 

AC e = mACt = 
2.66 X 10-* r>) faradays 

A schematic diagram of the com- 
puter circuit is shown in Fig. 1. The 
body of the computer, a conventional 
RC network made of non-inductive 
resistors and oil-filled condensers, is 
shown with each of the elements 
connected to the bars of a large 
commutator. The voltage source is a 
180 v "B" battery, chosen to give full 
scale deflection of the oscilloscope - 
beam when connected directly to the 
deflection plates. 

To simulate the movement of a 
heat source, the voltage source is 
connected to a rotating brush which 
contacts the commutator bars. A sec-' 
ond brush, rotating at a high speed 
relative to the primary voltage source, 
is used to transfer the condenser volt- 
age at each element to the oscillo- 
scope plates. By means of an adjust- 
able selector switch synchronized 
with the source rotor, a triggering 
signal is sent to a commutator bar 
adjacent to the last RC element to 
initiate the oscilloscope sweep. This 
results in a visual representation of 
the voltage profile in the computer 
network for the selected location of 
the simulated heat source. A photo- 
graph of the assembled computer and 
associated equipment is shown in 
Fig. 2. 
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To operate the computer, the rate 
of travel of the source rotor is se- 
lected by changing 'the gear ratio on 
the drive motor. After a suitable 
warm-up period, the positions of the 
trigger selector switch are adjusted 
and synchronization is checked with 
the oscilloscope sweep circuits. The 
camera is attached to the oscillo- 
scope and the network voltage pro- 
file photographed at the preselected 
positions. 

If the computer is used to simu- 
late a constant temperature source 
moving through a radial formation, 
a temperature profile similar to that 
shown in Fig. 3 is observed on the 
oscilloscope screen. The vertical dis- 
tance between grid lines represents a 
temperature difference of 100°F.The 
curvature of the grid lines calibrates 
the beam deflection for the entire 
screen. The voltage of each segment 
of the computer network is shown as 
a spot on the screen at a vertical 
position corresponding to its numeri- 
cal value. 

This analog computer can be 
adapted to the solution of heat trans- 
fer problems involving steady state 
or non-steady state conditions and 
moving or stationary heat sources. In 
order to check the network design, 
analog solutions were determined for 




Fig. 2 — Assembled Computer and 
Associated Equipment. 



problems involving a stationary heat 
source for which the theoretical 
curves could be obtained by conven- 
tional mathematical methods. In gen- 
eral, it was concluded that the net- 
work satisfactorily represented a cy- 
lindrical system and that the solution 
of non-steady state problems in a 
simulated infinite cylindrical medium 
was feasible for short times. For 
problems involving long computing 
times at a given position of the source, 
the error increases with time. 12 For 
such problems the accuracy could be 
improved by using a DC amplifier at 
each lump in the analog." ,14,15 How- 
ever, for problems involving a mov- 
ing source, the large additional ex- 
pense involved in such a refinement 
was not considered necessary. 

As described herein, the computer 
has been designed to study heat con- 
duction problems involving a mov- 
ing, constant temperature source in 
a radial formation. It should be 
pointed out, however, that the gen- 
eral design is applicable to a variety 
of heat transfer problems involving 
a moving source. With some modifi- 
cation' it would be possible to study 
non-steady state heat transfer by 
gases, liquids, and solids, either in- 
dividually or as an interacting sys- 
tem, in media of different geomet- 
rical shapes. In addition, the moving 
source rotor may be designed to rep- 
resent a heat source with any of a 
number of unique properties. 

ILLUSTRATIVE PROBLEMS 

To illustrate the type of problem 
which can be studied with this com- 
puter, several assumed modes of 
travel of the heat source will be dis- 
cussed. It should be emphasized that 
these problems are intended to be 
illustrative only and are not neces- 



sarily representative of a practical 
field operation. 

For the following illustrations, it 
is assumed that a constant tempera- 
ture heat source is produced in the 
formation. Although a number of 
possible mechanisms for the produc- . 
tion of such a heat source may be 
visualized, for convenience a mecha- 
nism similar to that postulated by 
Kuhn and Koch 4 has been assumed. 
According to their hypothesis, a com- 
bustion zone is produced by burning 
residual carbonaceous material de- 
posited out of the native crude oil 
ahead of the high temperature zone. 
In the idealized case, the tempera- 
ture of this combustion zone could 
be maintained at a constant optimum 
value by adjustment of the rate of 
air injection into the formation. 
Although reference will be made 
throughout the ensuing discussion 
to the "combustion front" and the 
associated fuel requirements, any 
means of obtaining a constant tem- 
perature front would result in identi- 
cal heat transfer characteristics for 
the particular conditions postulated. 

In one set of tests, a constant 
temperature source was moved radi- 
ally through the formation at various 
constant speeds which ranged from 
0.874 ft/ day to 0.097 ft/ day. Tem- 
perature profiles were photographed 
at convenient positions. In Figs. 4, 
5, 6, and 7, temperature data from 
the above tests have been plotted on 
a semi-log scale to show the com- 
puted temperature distributions for 
different positions of the heat source. 
It is immediately apparent that (1) 
at a constant rate of travel, preheat- 
ing of the formation ahead of the 
source isotherm increases with dis- 
tance of travel and (2) at a given ra- 
dius, preheating increases for slower 
rates of travel. 

These temperature distribution 
curves may be used to calculate the 
theoretically required amount of resi- 
dual fuel to maintain a self-propagat- 
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Fig. 3 — Typical Voltage Proftle 
in Analog Computer Network. 
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Fig. 4 — Temperature Distribu- 
tion Ahead of Combustion Front. 

Conduction only in radial for- 
mation. Rate of movement oT 
Trout : 0.874 ft/day. 

ing combustion front at a given rate 
of advance. Tn order to make such a 
calculation, it is necessary to know 
the instantaneous heat loss from the 
combustion front. If this heat loss 
can be supplied by burning the avail- 
able fuel, the temperature of the 
combustion front will remain con- 
stant and the source will be self-per- 
petuating as long as sufficient air is 
supplied. If there is insufficient fuel, 
ihe temperature will decrease and, 
under certain conditions, the process 
will, cease. 

The instantaneous heat loss can 
re determined from the equation, 
q = kA dT/dr. The term, dT/dr, is 
simDly (he instantaneous temperature 
gradient at the heat source and can 
be readily determined directly from 
the figures as 

dT/dr = 



Then the fuel required to maintain 
the combustion front will be 

k dT/dr 



W — - > w ' f V§ ib/cu ft of formation 
.v dr/dt 

where x is the heat of combustion of 
the residual fuel (assumed as 17,000 
Btu/lb) and dr/dt is the rate of 
movement of the front in feet per 
hour. 

Through use of the above method. 
Fig. 8 has been plotted to show the 
residual fuel required to maintain a 
moving heat source (temperature = 
I.O0O°F above ambient) in a forma- 
tion in which heat is transferred ra- 
dially by conduction alone. 

A suitable change in conversion 
factors (/ = 0.18, n = 0.006, m = 
6.08 X 10 '\ r - r ,/V2), can be 
made to show the effect of changing 
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Distance ahead of combustion front-ft. 

Fig. 5 — Temperature Distribu- 
tion Ahead of Combustion Front. 

Conduction only in radial for- 
mation. Rale of movement of 
front: 0.353 ft/day. 

the assumed value of the conduc- 
tivity. In Fig. 9 the fuel requirements 
to maintain a moving combustion 
front are indicated for k = 0.8, AT 
= 1,000. 

Similarly, through the choice of a 
different set of conversion factors 
(I = 0.30, n = 0.006, m = 6.08 X 
10" J \ r = r»/V2) f the effect of a 
change in the temperature of the 
combustion front can be found. Fig. 
10 shows the fuel requirements if 
k = 0.8, aT = 600. 

Tn general, Figs. 8, 9, and 10 show 
that the fuel requirements are higher 
for the slower rates of travel of the 
source. Inasmuch as the physical na- 
ture of the assumed process makes 
available only a given amount of 
fuel, which may be affected by for- 
mation and oil characteristics, there 
should be a minimum rate of ad- 
vance possible for a self-propagating 
combustion front. The curves indi- 
cate that this critical rate will de- 
crease with distance from the well- 
bore. This characteristic is a result 
of the greater preheating by conduc- 
tion as the combustion front travels 
away from the wellbore and should 
be true when the additional effect of 
iniected air is considered. 

From the data presented above, it 
is apparent that the differences be- 
tween the actual value of k and that 
used in the comDUter are not critical 
in determining the reauired fuel con- 
tent, provided k is of the right order 
of magnitude. Of greater importance 
in using the data is the value as- 
sumed for the temperature at the 
combustion front, inasmuch as the 
required fuel content was shown to 
be directly proportional to the tem- 
perature. This indicates one area in 
which experimental data from labo- 
ratory tests would be desirable. 



In order to maintain a constant 
rate of advance of the combustion 
front as shown in these figures, it 
would be necessary to provide for an 
increasing volume rate of injection 
air as the front travels away from 
the wellbore. To determine the re- 
quired air rates, it would be neces- . 
sary to know the position of the 
front at any time. This is not likely 
to be practicable for field operation. 
A simpler operational method would 
be to inject air at a constant rate. If 
both the fuel supply per unit volume 
of formation and the air injection 
rate are assumed to be constant, it 
would appear that because "of geo- 
metrical factors the minimum rate of 
movement of the combustion front 
should be inversely proportional to 
the radial position. 

It is possible to study the feasibil- 
ity of this mode of operation by 
either of two methods. In one case, a 
family of curves, each of which is 
proportional to 1/ radius, can be ob- , 
tained with a constant rate drive. In 
the second case, the computer can be 
set up so that the source rotor can be 
moved at a rate proportional to 
1 / radius. 

To illustrate the first method, the 
family of curves in Fig. 8 has been 
used to examine a mode of travel in 
which the rate of movement of the 
combustion front is equal to 8/ ra- 
dius in feet per day. Then for the 
constant rates of 0.874, 0.353, 0.211 
and 0.097 ft/day, the front would 
be located at approximately 9, 23, 38 
and 83 ft, respectively, in order to 
he traveling at the designated rate. 
If these points are plotted on Fig. 8 
they will all fall at a constant fuel 
requirement of approximately 2.3 
Ib/cu ft of formation, indicating that 
the front would be self-propagating 
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Fig. 6 — Temperature Distribu- 
tion Ahead of Combustion Front. 

Conduction only in radial for- 
mation. Rate of movement of 
front: 0.211 ft/dny. 
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under the above conditions. 

H the computer is set up so that 
temperature profiles are ohtained at 
radii, r if at the corresponding rate, 
8/n, the required residual fuel can 
be calculated from the computed 
temperature gradients in the manner 
described previously. In this case a 
single curve would* be obtained as 
shown in Fig. I I. Here asain the re- 
quired residual fuel is found to be 
2.3 Ib/cu ft of formation. 

in general, it may be concluded 
I hat for this mode of travel the nu- 
merical value of both the tempera- 
ture gradient and the radial velocity 
of the combustion front must be 
changing at the same rate. The results 
indicate also that in the radial sys- 
tem described a constant rate of fuel 
consumption will result in a constant 
temperature heat source only if the 
rate of movement is inversely pro- 
portional to the radial position. 

In the above discussion, heat trans- 
fer by the injected air has been neg- 
lected. It is possible, however, to use 
the available data to estimate the 
minimum fuel requirements in a sys- 
tem in which no heat is retained be- 
hind the combustion zone. This can 
be done by assuming that all of the 
heat stored in the rock as the mov- 
ing source reaches and passes a given 
radius will be transferred instantane- 
ously by the flowing air stream back 
to the combustion front. Thus the 
minimum fuel requirement in such a 
system can be approximated simply 
by subtracting the equivalent fuel 
content of the corresponding volume 
increment of heated rock from the 
residual fuel content calculated from 
the data shown in Figs. 4, 5, 6, and 
7. This has been done for a radial 
thermal system in which k- 1.6 and 
AT = 1,000 and the calculated 
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Fig. 7 — Temperature Distribu- 
tion Ahead of Combustion Front 

Conduction only in radial for- 
liuititm. Rale of movement of 
front: 0.097 ft/diiy. 
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RADIAL DISTANCE FROM WELL 90RE -FEET 

Fig. 8 — Fuel Required to Main- 
tain a Moving Combustion Front. 

Thermal conductivity: 1.6 Btu/ 
ft r hour ^F/fu Combustion front 
1,000 C F above Ambient Tem- 
pera lure. 

curves have been plotted in Fig. 12. 
The values for the residual fuel con- 
tent shown in Figs. 8 and 12 repre- 
sent two extremes of operation. In 
the one case (Fig. 8), it is assumed 
that all of the heat stored in the rock 
behind the combustion front remains 
unavailable for raising the tempera- 
ture of the formation. In the second 
case (Fig. 12), all of this heat is uti- 
lized to the fullest extent. 

As a final illustration of the use- 
fulness of the analog, there is shown 
in Fig. 13 the analog computation 
of the temperature decay in a radial 
formation after a heat source has 
stopped moving. In this problem, the 
combustion front was moved through 
the formation at a constant rate of 
0.353 ft/day to a radius of 40 ft and 
then stopped. This might correspond 
to a system (with heat flow by con- 
duction only) in which air injection 
has stopped because of equipment 
failure. From the slow decay of tem- 
perature shown, it might be con- 
cluded that once a hot zone has been 
well established, re-ignition of the oil 
sand following a prolonged shutdown 
should not be a serious problem. 

This problem is similar in nature 
also to that investigated mathemati- 
cally by Jenkins and Aronofsky. 10 The 
temperature increase after 840 days 
of 70° F at a radius of 100 ft com- 
pares reasonably well with their theo- 
retical value of 75° F after two years. 

From the above illustrations, it 
would appear that the computed tem- 
perature distribution curves should be 
useful for studying various factors, 
such as the best start-up technique 
and the optimum mode of travel of 
the heat front. It is obvious, however, 
that the simplified mechanism set up 
for the foregoing problems does not 
simulate closely the heat transfer 
process existing in a practical field 
operation. In order to obtain the 
maximum utility from this machine, 
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Fig. 9 — Fuel Required to Main- 
tain a Moving Combustion Front. 

Thermal conductivity: 0.8 Biu/ 
fr hour °F/ft. Combustion front 
1,000 °F above Ambient Tem- 
perature. 

it would be desirable to modify it to 
include additional factors such as 
heat transfer by the injected air and 
by the movement of hot fluids ahead 
of the combustion front. It seems ap- 
parent that modification of the com- 
puter to study heat transfer by vari- 
ous assumed mechanisms must be 
coordinated closely with results ob- 
tained from laboratory and field ex- 
perimentation. 

CONCLUSIONS 

A design has been presented for 
an electrical analog computer which 
can solve non-steady state heat con- 
duction problems pertaining to an ex- 
tensive radial system containing a 
moving, cylindrical heat source. 

It has been the intent of this paper 
to show the type of problem the ana- 
log computer can be adapted to 
study, and to indicate how its use 
can be integrated with the results 
obtained from laboratory data to in- 
crease our understanding of possible 
field applications of oil recovery by 
thermal methods. A method has been 
shown for using the temperature dis- 
tribution curves obtained with the 
computer to calculate the fuel re- 
quired to maintain a self-propagating 
heat source. 

Inasmuch as the computer design 
presented herein takes into account 
only those heat losses resulting from 
radial conduction, the data should 
not be construed as applying directly 
to field operations. Modification of 
the computer to include other heat 
transfer mechanisms depends pri- 
marily upon an understanding of 
the particular thermal oil recovery 
method being studied. Laboratory 
scale studies of. the thermal process 
itself aid in determining the most 
satisfactory mechanism to use. 

Because of the complex nature of 
thermal oil recovery processes, it 
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Fig. 10— Fuel Required to Main- 
tain a Moving Combustion Front. 

Thermal conductivity: 0.8 Bin/ 
ft 2 hour °F/ft. Combustion front 
600 o F nhovc Ambient Tempera- 
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Fig. 11 — Fuel Required to Main- 
tain a Moving Combustion Front. 

Thermal conductivity: 1.6 Btu/ 
ft 2 hour °F/ft. Combusion tem- 
perature 1,000 °F above Ambient. 



would appear that continued efforts 
of research workers in many differ- 
ent groups will be required to eval- 
uate the numerous possible ramifica- 
tions in a reasonable time. It is hoped 
that this paper will stimulate the pub- 
lication of other work on thermal 
recovery processes. 
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NOMENCLATURE 

Thermal units 

i? t = thermal resistance, 

°F hour/Btu 
C t = thermal capacity, Btu/°F 
Qt =heat stored in C„ Btu 
q — hourly heat flow, Btu/hour 
A T = rise in temperature of C t 
due to Q t , °F 

combustion front moves 0353 ft/ day 
for 40 ft; then stops 




radial distance in feet 

Fig. 13 — Conduction of Heat 
After a Combustion Front 
Has Been Stopped. 



r t = distance in the thermal cir- 
cuit (radial), ft 
,A n = average cross sectional area 
of each increment Ar„, 
sqft 

N — number of the element 
(Ar) under considera- 
tion 

r t = time in the thermal circuit, 
hour 

k 

oc = — = thermal diffusivity, 

CD 

sq ft/ hour 
c — heat capacity per unit 

mass, Btu/lb °F 
k = thermal conductivity, Btu/ 

sq ft hour °F/ft 
p = density of formation, lb/ 
cu ft 
Electrical units 

Re = electrical resistance, ohms 
C e = electrical capacity, farads 
Q a — charge stored in C B , cou- 
lombs 

/ = electrical current, amperes 
&V — rise in voltage of C 0 due 
to 2-» volts 

/„ = time in the electrical cir- 
cuit, seconds 

r c = distance in the electrical 
circuit (radial) , ft 
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The authors are to be compli- 
mented for their ingenious solution 
to a problem having an important 
bearing on future development of 
thermal oil recovery processes. The 
interpretation of results from this 
heat transfer study relative to an in- 
situ combustion oil recovery process 
illustrates the utility of this type of 
theoretical work. Undoubtedly, their 
electric analog computer will have 
many other uses in reservoir tech- 
nology. 

With regard to the results of the 
study, it should be emphasized that 
for the more conventional well pat- 
terns, radial flow conditions exist for 
only 20 to 30 per cent of the distance 
from an injection well toward the 
nearest production well. The results 
from this study are, therefore, most 
applicable and important in consid- 
eration of the ignition phase of this 
type of thermal oil recovery process. 
As the authors implied, the fuel re- 
quirements presented are based on 
assumed thermal properties and an 
idealized mechanism. Actual fuel con- 
sumption can differ greatly from the 
theoretical values obtained from this 
study. Nevertheless, the trends illus- 
trated by the authors permit impor- 
tant observations. 

For the authors 1 case of a constant 
combustion front velocity, the mini- 
mum fuel requirement appears espe- 
cially critical when the burning front 
is near the injection well and/or 
when the front is moving at a slow 
rate. In contrast, the results for the 
case where the velocity varies in- 
versely with radius of the burning 
front indicate a fuel requirement that 
does not vary with position of the 
burning front. While this result tends 
to contradict the conclusions drawn 
from the case of constant velocity, 
it is believed that most actual opera- 
tions will fall between these condi- 
tions, and the fuel requirement will 
he substantially higher near the in- 
jection well. 
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During investigation of this ther- 
mal recovery process, we computed 
fuel requirements by an analytical 
solution of substantially the same 
heat conduction problem. Two cases 
were considered: (1) heat conduction 
from an infinite cylindrical source 
moving radially in an infinite medium 
of constant thermal properties, and 
(2) the same problem with the in- 
clusion of the boundary conditions 
resulting from consideration of the 
ignition wellbore. 

The solutions to Case (1) for a 
burning front velocity inversely pro- 
portional to the radial position of the 
front showed the fuel requirement to 
be entirely independent of position 
of the front. And for the same ther- 
mal conditions assumed by Vogel and 
Krueger, computed results were in 
good agreement with the values from 
the electrical analogy. A fuel require- 
ment of 2.4 lbs/cu ft was obtained 
as compared with 2.3 reported by 
Vogel and Krueger. This establishes 
a direct verification of the electrical 
analog for a moving source. 

A solution for Case (2) at a radial 
position of 0.5 ft indicated a slightly 
higher fuel requirement of 2.7 lbs/cu 
ft. But this is attributed to the effect 
of the wellbore. 

The analytical solutions mentioned 
above are as follows. 

1. Case (1) — Constant velocity 
of source (r' = Ut') . 

T - T n — 



e A er (t - £') 



U_q( f 
« 

r ( rUt' ) 

M^cT^T)}^- •••(!> 

Where T 0 is the initial temperature, 
T is the temperature at radius r at 
time t, q is the rate of heat genera- 
tion per square foot of frontal area, 
U is the constant velocity of the 



burning front, k is the thermal con- 
ductivity, t is time, and t' is the time 
the burning front is located at radius 
r', oc is the thermal diffusivity, and 
/ 0 is Bessel's function of the first kind 
of an imaginary argument. If Equa- 
tion 1 is integrated for r = r' at f', 
the fuel requirements may be found 
from the expression: 

lb fuel _ q 

Th' ' ' ' 



(2) 



cubic foot 

where s is the heat of combustion of 
the fuel. 

2. Case (1) — Velocity of source 
is inversely proportional to the radius 
Or" = 2Vt). 



T - T a = 



qV2V( y/f 



2k 



V V* 



(r» + 2VV) 
4 oc (t - f) 



/ V y/2Vf \ 

M 2cc(<-o r' • ■ ■ • (3) 

where V = Ur' = constant. . (4) 
The fuel requirements may be found 
by integrating Equation 3 for r = r' 
at t' and using the expression: 
lb fuel ^ qr f 
cubic foot ~~ 7F * * ' * ^ 

3. Case (2) — Wellbore effect in- 
cluded = 2 

cpV 

i 



(£)■ ■ • 



xdxd 
where 
C 



(6) 



... (7) 

In Equations 6 and 7, c is the spe- 
cific heat of the saturated sand, p is 
bulk density, r w is the radius of the 
wellbore, /„ and Y„ are Bessel func- 
tions of the first and second kind, 
respectively, of order zero, and x is 
a variable of integration. 

A further refinement of the prob- 



lem treated by authors Vogei and 
Krueger is the consideration of verti- 
cal heat losses. This would require 
the expensive modification of rebuild- 
ing the RC network in two dimen- 
sions, and considering a burning front 
of finite vertical dimension. Vertical 
heat loss can be important to the 
maintenance of self-sustained com- 
bustion when the sand section is thin 
and the well spacing is large. If ver- 
tical heat losses were considered, it 
would also be desirable to account 
for the lower thermal diffusivity of 



the "burned" sand section of the 
reservoir. 

We agree with the authors' state- 
ment that the data developed by 
them in this paper should not be 
construed as applying direcdy to field 
operations, inasmuch as the front is 
known to progress over a range of 
temperatures, velocities, and fuel con- 
sumptions in accordance with incom- 
pletely defined rate phenomena oc- 
curring in a reservoir. Many chal- 
lenging problems along these lines 
remain to be solved. 
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We should like to thank Rodman 
Jenkins and H. J. Ramey, Jr., for 
their pertinent comments. It is grat- 
ifying to see the good agreement re- 
ported between the analytical and 
the analog solutions to the prob- 
lem of determining the fuel require- 
ment in a radial system in which 
the burning front velocity is in- 
versely proportional to the radial 
position. This same agreement has 
been noted by R. G. Hawthorne of 
our laboratory. 

The presentation of the equations 
used in studying the simplified prob- 
lems discussed in the paper serves 
to illustrate the somewhat compli- 
cated nature of analytical solutions 
to heat transfer during thermal oil 
recovery processes. Equations 3, 6, 
and 7 illustrate the complexity added 
to the simple case by considering 
the wellhore effects. In the com- 
puter, wellbore effects are au- 
tomatically taken care of in the 
design of the electrical lumps near 
the wellbore. It would be interesting 
to see the analytical solutions to the 
more complex heat transfer problem 
involving the additional factors of 
vertical heat losses and heat trans- 
fer by the driving gases and con- 
densing fluids. 



Tn general, analytical solutions of 
this nature are tedious and time con- 
suming unless attacked with modern 
digital computer techniques. In any 
case, considerable expense is in- 
volved in arriving at numerical an- 
swers. It would be of considerable 
value to the industry to have further 
work published on solutions to the 
more complex problems. On the 
basis of our present thinking, we 
feel that most such problems can be 
treated somewhat more simply and 
inexpensively by the modification of 
the analog computer. 

Jenkins and Ramey have sug- 
gested that vertical heat losses be 
considered as a refinement to the 
problems treated in the paper. As 
they have noted, vertical heat losses 
can be important when considering 
the problem of maintaining a com- 
bustion front in thin oil zones. We 
believe, however, that the applica- 
tions in which vertical heat losses 
are the controlling factor are likely 
to be the exception rather than the 
rule. For thicker zones there is some 
evidence from analytical studies that 
the effect of vertical heat transfer 
on fuel requirements is considerably 
less than the effect of heat trans- 
fer by injected air (or inert gases) 



and possibly by the flowing liquids. 
The effects of both of these latter 
factors can be studied rather easily 
by inexpensive modifications of the 
analog computer. 

As Jenkins and Ramey have em- 
phasized, radial heat flow conditions 
exist for only 20 to 30 per cent of 
the distance from the injection well 
toward the producing well. At the 
moment, we believe this region to 
require the most critical considera- 
tion when attempting to evaluate the 
feasibility of starting a successful 
thermal oil recovery process. How- 
ever, we feel that a study of heat 
transfer during later stages of the* 
process can be achieved through 
relatively simple modifications of 
the computer. It appears that ana- 
lytical solutions involving interfer- 
ence between wells might become 
quite complex. 

It is obvious from a consideration 
of both the comments of Jenkins 
and Ramey and the authors* reply 
that there is much to be studied and . 
learned in this field. Both analytical 
and analog studies can be helpful in 
solving the many problems involved. 
However, as mentioned in the paper, 
their maximum utilization will be 
achieved through correlation witi., 
laboratory and field studies. *** 
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